Spin-orbital coupling effect on Josephson current through a superconductor 
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We study spin-orbital coupling effect on the Josephson current through a superconductor (SC) 
heterojunction, consisting of two s-wave superconductors and a two-dimensional electron gas (2DEG) 
layer between them. The Rashba-type (RSOC) and/or Dresselhaus-type (DSOC) of spin-orbital 
coupling are considered in the 2DEG region. By using the lattice Bogoliubov-de Gennes equation 
and the Keldysh formalism, we calculate the DC supercurrent flowing through the junction and 
find that the critical current 7c exhibits a damped oscillation with both the strength of SOC and 
5IIJ , the layer length of 2DEG; especially, the strength ratio between RSOC and DSOC can also induce 

' switching between the state and the tt state of the SC/2DEG/SC junction as well. This O-tt 

transition results from the fact that SOC in a two-dimension system can lead to a pseudo-magnetic 
00 . effect on the flowing electrons like the effect of a ferromagnet, since the time reversal symmetry of 

the system has already been broken by two SC leads with different macroscopic phases. 

o 
o 

' The issue of Josephson current sign reversal in a superconductor/ferromaget/superconductor (SC/FM/SC) het- 
erojunction has recently drawn a great deal of attention[l|, [3, H due to its experimental observation and potential 
^ ■ applications to spintronics and quantum computing. SIhjHj Q The physical origin of the Josephson effect is the 
^ • ' breakdown of time reversal symmetry (TRS) in SC/SC junctions due to the SC macroscopic phase difference. Between 
. the two SC leads in the junction is either an insulator or a "weak link" such as normal metal (NM), semiconduc- 
^ ' tor, and so on. The DC Josephson effect can be understood by the Andreev reflection processes of quasiparticlesQ 
I . with energy smaller than the superconducting energy gap, an electron impinging on one of the interfaces is Andreev 
■ O ' reflected and converted into a hole moving in the opposite direction, thus generating a Cooper pair in an SC; this 
C , hole is consequently Andreev reflected at the second interface and is converted back to an electron, leading to the 
destruction of a Cooper pair in the other SC. As a result of this cycle, a pair of correlated electrons are transferred 
from one SC to another, creating a supercurrent flow across the junction. [l| 

When an FM is inserted between the two SC leads, the current-carrying Andreev bound states are spin-splitted 
and shifted in an oscillatory way. FM favors to align the electron spin due to the FM exchange interaction, whereas 
the ordinary spin singlet Cooper pair consists of two spin-antiparallel electrons; thus FM tends to destroy the su- 
perconductivity and the Josephson current in SC/FM/SC junctions would be depressed compared with SC/NM/SC 
junctions. Furthermore, the Cooper pair is associated with a nonzero momentum due to spin splitting from the FM 
exchange interaction and the Josephson current exhibits a damped oscillating behavior with the FM layer length or 
the FM exchange strength, i.e. , the O-tt transition. [H [HI, [H, [ll [H, [H, [M E H [H 113, ISl, [H HI The 7r-state of 
■ the SC/FM/SC junction is the Josephson current flowing in the direction opposite to the phase difference between the 
OQ ■ two SCs. To realize systematically the O-tt state switching in a single SC/FM/SC junction, it needs to change either 
. the length of the FM layer or the FM exchange strength; however, as a matter of fact, it is very difficult to do this in 
^ ' an experiment. Since the tunability of a system is of importance for experimental observation, some alternatives have 
, been proposed; for example, the SC/FM/FM/SC junction with noncoUinear magnetizations has been investigated by 
' Pajovic et al.W^ and the O-tt transition can be found by modulating the relative direction of the FM moments. Dolcini 
JH and Giazotto[25| proposed to use Aharonov-Bohm interferometry to realize a fully controllable Josephson junction. 
■ - - The Josephson current flowing through a two-dimensional electron gas (2DEG) region with Rashba (RSOC) and/or 

Dresselhaus (DSOC) spin orbital coupling has also been studied, because the pseudomagnetic field from SOC may 
lead to the same effect on the Josephson current as FM. j26i. i27i. |28. 129] Several studies[26,, ,27] have shown that in the 
one-dimension case the SOC can hardly exert any effect unless the Zeeman splitting from an external magnetic field 
is included. It was generally argued that SOC does not break the TRS of the system and cannot function like an FM. 
However, Dell' Anna et aL[2^ argued that TRS has already been broken by the supercurrent and found that SOC can 
have a huge effect on the Josephson current in a quantum dot system, although they did not find the O-tt transition 
behavior since only two electron levels are considered in the quantum dot. 

In this work, we restudy the SOC effect on the DC Josephson current through a SC/2DEG/SC junction in the 
clean limit and with the multilevels in the 2DEG being taken into account. When TRS of the system is broken by 
the two SCs with different macroscopic phases, the pseudomagnetic field from SOC should cause an effect (hereafter 
referred to as pseudo-magnetic effect) on the electron transport property, resembling the exchange field effect in a 
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FM; however, for ID transport, the pseudo- magnetic effect from SOC should disappear. Based on this analysis, we 
calculate the Josephson current flowing in a SC/2DEG/SC junction by using the discrete BdG equation and Keldysh 
Green's functions. It is shown that the Josephson current exhibits a damped oscillation as a function of the strength 
of SOC, the length of the 2DEG layer, and the strength ratio between RSOC/DSOC. This O-tt transition is the same 
as that found in SC/FM/SC junctions. However, the O-tt transition induced by SOC is much easier to be realized in 
experiment since the strength of SOC can be modulated to a large extent by an external electric field perpendicular 
to the plane of 2DEG.[3H 

We consider a clean SC/2DEG/SC heterojunction, consisting of two semi- infinite SC leads and a 2DEG layer with 
a length d between them, as shown in Fig. 1. The interfaces between 2DEG and SCs are set at a; = and x = d and 
the current flows along the x-direction. In the 2DEG layer both RSOC and DSOC are considered, which come from 
the structure asymmetry of the materialjH, [s^ and are expressed as 



a 13 

T^SOC = -^{f^yPx - CTxPy) + -^{(^xPx " O-yPy), (1) 

where a and f3 are the coupling constants of RSOC and DSOC, repectively, and RSOC can be directly modulated by 
an external electric field perpendicular to the 2DEG plane, |3lJ and Py are the two components of the momentum 
operator p, ax and ay are the Pauli matrices. The difference between these two SOCs is the different directions of 
the pseudomagnetic fields from SOC for the same electron momentum. By introducing a ratio angle 6, the equation 
above can be rewritten as 

'HsOC = -^{^yPx - ^xPy)- (2) 

/ e-''^ \ f -ie^^ \ 

Here j — + /J^, a = ■^cosO, and (3 = j sin 9, a^ — I ^ 1 , and ay = i j^^-ie q ] ■ ^he angle 9 denotes 

not only the strength ratio between RSOC and DSOC, but also the pseudomagnetic field direction for the mixture of 
RSOC and DSOC, which is similar to the rotation angle of the spin axis.[30t as can be seen in the renormalized ax 
and ay. 

The two SC leads are described in the framework of standard BCS formalism with a macroscopic phase difference 
(p. In the absence of external applied field, the mean-field Hamiltonian of the SC/2DEG/SC junction reads [33| 

Ti. — Hl + Ti-R + Ti2DEG + Wt, (3) 

^L(fl) = '^{^tja- - tJ')Cjj„Cija - t'^iCj+lja'^tj'y + Cjj+l^a'^ija + C.C.) 

ija 

-^(Ai(^)4^4^ + A*(«)C,,xa,T), (4) 



'H2DEG — '^{^ija ~ n)C}j^Cija- ~ t^^iCj+ijaCijo- + Cjj_^_^ ,^Cija + C.c} 
ija ijcr 



-tso {G}+ij„{iay)cra'Cijcr' — Clj^-^^^iiax) aa' Cijai -(-C.C.}, (5) 



nr = J2(*'lCI^lC,.,2Deg + t'jiClj,C)a,2DEG + c.c). (6) 

JO" 

Here Hl and Hr are respectively the discretized BdG equation of the left and right SC leads; Ti.2DEG is the 
Hamiltonian of the 2DEG region, i.e., a free electron model with SOC included; Ht is the tunneling Hamiltonian 
between the 2DEG and the left and right SC leads. Cjj^{Cijcr) is the creation (annihilation) operator of an electron at 
site (ij) with spin a =ti, Siju = 4t is the site-energy, t = /2'ma^ is the hopping energy with the lattice constant a 
, and tso — 7/2a denotes the SOC strength in the lattice representation. Ai^(^r-) = | Ai(j^) je^***/^ is the pair potential 
in the left (right) SC lead with a macroscopic phase ±0/2, fi is the chemical potential which is constant everywhere 
in our model due to zero applied bias in the junction, ^'^(j^) is the hopping strength between the left (right) lead 
and 2DEG, which is independent of spin so that no spin flip effect occurs when quasiparticles tunnel through the 
interfaces; it can also represent the strength of the interface barrier between the SC lead and the 2DEG region. 
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Although we consider a clean system without impurity, the Hamiltonian above is not limited to this system and 
can describe systems with other band structures as well as different pair potential symmetry. Along the y direction of 
the junction in Fig. 1, the translational symmetry is preserved so that the electron momentum ky is a good quantum 
number and summation in Eq. 3 over j sites can be transformed to the /c-space, e.g., 

'H2DEG = ^(4i - - 2tcoskya)ClC - + c.c) - 

^^{tsoC}^-i^„{iay)Cia' + c.c. + 2tso sin kyaC}„{ax)Cia'). (7) 

Other components of the Hamiltonian (Eq. 3) can be treated in the same way, thus the numerical calculation will be 
reduced greatly compared with a finite size (along the y-direction) system. We focus on the supercurrent through the 
interface between the SC and 2DEG, thus, the current density operator is given by = e^^^ = jj-[Nl, Ht] where 
Nl is the electron operator in the left SC, and after commutation the steady current density reads [35<] 

Jl = TTTp I ^{*'LGt,2DEGi^^ h) ^ (4)*G'^D£;G,l(^' fcy)}ll+33- (8) 

G'^{E,ky) is the lesser Green's function defined in 4 x 4 Nambu0spin space, Gf^ j^,{t,t') = i{Cj^,{t')Cia{t)) with 
(. . .) denoting the quantum statistical average. W is the transverse width of the junction, the subscripts 11 + 33 
stands for the summation over the spin-up and spin-down components of the electron current, i.e., the (l,!) and 
(3,3) matrix elements. Since we consider the DC Josephson effect and no bias is applied on the system, the lesser 
Green's function could be obtained by the equilibrium equation G^{E) — {G°-{E) — G''{E))f{E) with f{E) being 
the Fermi-Dirac distribution function, G""^"-* is the retarded (advanced) Green's function, which fulfills the Dyson 
equation, e.g., 

where g^*^"'' is the surface Green's function of the left SC lead and Gj^^q is the coupled Green's function of the 
2DEG, Htl is the left tunneling matrix of Eq. (3d) in the Nambu{^spin space. The retarded Green's function can be 
calculated by direct matrix inversion G2DEG = ^ H2DEG ~ ~ '^r)~^ the recursive method, [36] where H2DEG 
is the 2DEG Hamiltonian Eq. (3b) in the Nambu{^spin space and S^^-, is the self-energy from the left (right) SC 
lead. 

In the numerical calculation, the hopping energy t is set as the energy unit t — I, the chemical potential takes 
fi — 2.0t and the pair potentials in the two SC leads are assumed to be equal, lA/^f/jjl = O.Olfi. The Josephson current 
is calculated at zero temperature based on Eq. 5 and the current phase relations are shown for different strength 
ratio angles 9 in Fig. 2. The hopping term stands for the interface transparency between the SC leads and 

2DEG; t']^ — ^'r ~ t is set in the calculation to denote full transparency of the interfaces so that deviates from the 
sinusoidal relation as shown in Fig. (2a). The 6 = n/A case (dashed-line) exhibits the coexistence of the and tt states, 
which suggests a O-tt transition of the junction should occur upon changing the angle 9 and other parameters. In 
Fig. (2b), it is expected that /(</>) can approximate a harmonic function for low interface transparency t'^ = t'j^ = 0.5t. 

In Fig 3, the critical Josephson current Ic is presented as a function of the angle 9 for different SOC strength tso and 
2DEG layer length d. The dips in these curves correspond to the O-tt transitions and the oscillation period decreases 
with the SOI strength. Changing 9 is equivalent to modulating the relative pseudomagnetic field direction, which is 
perpendicular to the electron-momentum direction for RSOC while for DSOC, the pseudomagnctization direction is 
complementary to RSOC by exchanging the spin axis.^SflJ Here the mixture of RSOC and DSOC is similar to the 
SC/FM/FM/SC junction with noncoUinear FM moments studied by Pajovic et aZ.,[23] who found Ic can oscillate 
with the angle between the two magnetizations. In Fig. 4 and Fig. 5, the critical current Ic is also shown to exhibit 
the O-TT transition by varying either the SOC strength tso (Fig. 4) or the 2DEG layer length (Fig. 5). The damped 
oscillation of Ic suggests that the 7r-state of the junction stems from the psuedo-magnetic effect of SOC, which is 
similar to an SC/FM/SC junction in terms of physical origin. 

The TT state of the SC/FM/SC junction results from the Cooper pair (coherent electron and hole in FM) being 
associated with a nonzero momentum due to the FM exchange interaction, which can lead to an Ic oscillating with 
the FM exchange strength as well as the FM layer length. For the studied SC/2DEG/SC junction, SOC maintains 
the TRS of the system itself on the whole because each electron sees a pseudomagnetic field B = ^(z x p) from Eq. 2; 
the system does not have any pseudo-magnetic effect, which resembles breaking the TRS, when a summation is taken 
over all electrons. However, in the present case the TRS of system has already been broken by the current, and the 
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SOC should have effects on the transport properties. By summing over those electrons (momenta) contributing to 
the current, the pseudo-magnetic effect of SOC appears; for example, the anomalous tunneling magnetoresistance 
(ATMR) observed in experiments [s^] is generally attributed to the SOC effect which is a pseudo-magnetic effect on 
the TMR. When a Josephson current is flowing through the SC/2DEG/SC junction, SOC can affect it definitely and 
Ic can exhibit a damped oscillation like in an SC/FM/SC junction. 

For one-dimensional SC/2DEG/SC junctions with SOC or only perpendicular injection considered [ky = 0), one 
cannot obtain O-tt transition as those shown in Fig. (3-5). The SOC cannot lead to a pseudo- magnetic effect in the ID 
case, because in this case the density of states is still spin-degenerate, unlike in the 2D case. [30] Therefore, ky ^ is 
a prerequisite for SOC exhibiting magnetic effect on electron transport and the summation over all transverse modes 
ky does not mean the effect of SOC would disappear, for the current is not simply an odd function of ky, at least 
in the clean limit. This is different from the semiclassical average over ky in Ref. 28. Compared with SC/FM/SC 
jucntions, the SOC strength as well the ratio angle are relatively easy to modulate by an perpendicular electric field, 
this feature may be of importance for experimental applications. For instance, Nitta et al. found the SOC constant 
can be modulated up to 50% in InAs material. [sH] 

In summary, we have studied SOC effect on the Josephson current through a SC/2DEG/SC junction based on 
the lattice BdG equation and Keldysh formalism. It was found that the critical Josephson current Ic can exhibit a 
damped oscillation with both the strength of SOC, the layer length of the 2DEG, and the strength ratio angle 9. The 
O-TT transition, same as that found in SC/FM/SC junctions, results from SOC induced pseudo-magnetic effect in a 
2D system when TRS of the system is broken by the supercurrent flow; but for ID junction the tt state is absent. 
The findings here can be observed experimentally by modulating the strength of SOC with external electric field. 
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Figure 1 Yang et al. 
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FIG. 1: The schematic of a two-terminal device in which two SC leads connect with the 2DEG, the shadow parts are the two 
interface barriers. Josephson current resulted from the macroscopic phase difference of the two SCs flows through the 2DEG 
along the x-direction. 




FIG. 2: Current phase relation I{<j}) for different barrier strengths, (a)i^ — t'^ — t and (b) t —'j^= t'n = 0.5t. Other parameters 
are d = 20a, tso = 0.05, = Q (sohd Une), 6 = 0.257r (dash Une), and 9 = O.Stt (dot-dash hne). 




FIG. 3: The critical current as a function of angle 6 for different SOC coupling strengths, tso = 0.05 (dash line) and t, 
(solid line), d = 20a. 



